Theory for superfluidity in a Bose system 
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We present a microscopic theory for superfluidity in an interacting many-particle Bose system 
(such as liquid *He). We show that, similar to superconductivity in superconductors, superfluidity 
in a Bose system arises from pairing of particles of opposite momenta. We show the existence of an 
energy gap in single-particle excitation spectrum in the superfluid state and the existence of a specific 
heat jump at the superfluid transition. We derive an expression for superfluid particle density Us as 
a function of temperature T and superfluid velocity Vs . We show that superfluid-state free energy 
density F is an increasing function of Vs (i.e., dF/dvs > 0), which indicates that a superfluid has 
a tendency to remain motionless (this result qualitatively explains the Hess-Fairbank effect, which 
is analogous to the Meissner effect in superconductors) . We further speculate the existence of the 
equation j — — AV x cj, where j — naVs is the superfluid current density, uj = V x Vs the superfluid 
vorticity, and A a positive constant (with the help of this equation, the Hess-Fairbank effect can be 
quantitatively described). 

PACS numbers: 67.25.D- 



I. INTRODUCTION 

We present in this paper a microscopic theory for su- 
perfluidity in an interacting many-particle Bose system 
such as liquid "^He. The theory is based on an assump- 
tion that particles of opposite momenta are paired in the 
superfluid state, and thus, is similar in many respects to 
the BCS theory of superconductivityji 

It is well known that there is a marked similarity be- 
tween liquid ''He II (the superfluid phase of liquid ''He) 
and superconductors, both being chiefly characterized by 
their ability to sustain flows of particles at a constant 
velocity without a driving force. ^'■^ However, unlike su- 
perconductors, for which there exists a successful micro- 
scopic theory, i.e., the BCS theory of superconductivity,^ 
a satisfactory microscopic theory for liquid *He II is 
still lacking, despite many efforts (for example, Refs. 

iiiES). 

Fundamental to the BCS theory of superconductivity 
is an assumption that electrons of opposite momenta and 
spins are paired in the superconducting state.— This as- 
sumption allows microscopic derivation of all essential 
properties of the superconducting state, such as the exis- 
tence of an energy gap in electronic excitation spectrum, 
a second-order phase transition (manifested by a specific 
heat jump at the superconducting transition), the Meiss- 
ner effect, and the Josephson effect. 

In this paper we show that it is also the pairing of 
particles of opposite momenta that is responsible for su- 
perfluidity in a Bose system. Namely, the cause for su- 
perconductivity in superconductors and superfluidity in 
liquid ^He II is indeed essentially the same, irrespective 
of the nature of the particles involved. 

Some previous attempts to develop a microscopic the- 
ory for superfluidity in liquid ''He II failed at the very 
start by assuming that the ground state of liquid He II 
is a Bose-Einstein condensate (for example, Ref. H)- As 
we will see in this paper, the ground state of a superfluid 



is not a Bose-Einstein condensate, but a state in which 
particles of opposite momenta are paired, similar to that 
of superconductors. 

Pairing of particles in a Bose system has been stud- 
ied by a number of authors (for example, Refs. [ol llOlllll . 
|l3 ). However, the authors did not treat properly self- 
consistency associated with pairing approximation, and 
thus, failed to establish a connection between pairing and 
superfluidity. 

In Sec. [ni we present the theory for the case where a 
superfluid is at rest, and show the existence of an energy 
gap in single-particle excitation spectrum in the super- 
fluid state, and the existence of a speciflc heat jump at 
the superfluid transition. In Sec. IIIIl we present the 
theory for the case where a superfluid current is present. 
We derive an expression for the superfluid particle den- 
sity as a function of temperature and superfluid velocity. 
We show that the superfluid-state free energy density is 
an increasing function of superfluid velocity, which indi- 
cates that a superfluid has a tendency to remain motion- 
less. This result provides a qualitative explanation for 
the Hess-Fairbank experiments^ in which a reduction of 
moment of inertia was observed when a rotating cylinder 
of liquid ^He was cooled through the superfluid transi- 
tion (this phenomenon, known in the literature as the 
Hess-Fairbank effect, is analogous to the Meissner effect 
in superconductors). We further consider how the Hess- 
Fairbank effect can be quantitatively described. A brief 
summary is given in Sec. IIVI 



II. SUPERFLUID TRANSITION 

We consider an interacting many-particle Bose system. 
We assume in this section that superfluid velocity Vs = 
(we will consider the case where Vg ^ in the next 
section). 

Similar to the pairing Hamiltonian in the BCS theory 



of superconductivityji we write the Hamiltonian of the 
interacting many-particle Bose system as 



the single-particle excitation energy in the superfluid 
state; and 



H = J2(^^ 



-m) 



t 1 V- 



2X!^''k'^k«-k«-k'ak', (1) 

kk' 



where ek is the normal-state single-particle energy, fi the 
chemical potential, Vkk' the pairing interaction matrix 
element, and aj^ and ak are Bose operators for a single- 
particle state of wave-vector k in the normal state and 
satisfy the commutation rule [ak,a]^,] = 5k, k' • 

This Hamiltonian can be diagonalized in essentially the 
same manner as in the BCS theoryiiiiiii^ Namely, we 
assume 



(a-kflk) 7^ 



(2) 



"-k 



^Ek/ksT __ ^ 



(9) 



the Bose function (the number of single-particle excita- 
tions of wave- vector k) . 

Coefficients Wk and Wk are found to satisfy the following 
relations: 



i-i'4 1 + ' 



kk =7^ T^-1 



2 V Sk 



(10) 



(11) 



in the superfluid state for a pair of (k) and (— k) par- 
ticles (where the angle brackets (• • •) denote a thermal 
average); treat (a_kak — (a-kOk)) as a small quantity 
so that terms bilinear in (o-kOk — (a-kflk)) can be ne- 
glected; define an energy gap parameter 



E^- 



k'(a_k'ak' 



(3) 



(because of the similarity between the present theory and 
the BCS theory of superconductivity, we will similarly 
refer to the quantity Ak as an "energy gap parameter" 
in this paper, although, as we will see below, it does not 
directly relate to an "energy gap" in the present theory) ; 
and apply a canonical transformationa l "'^'^'^^ 



and 



AkUkWk 



2^k 



(12) 



After the diagonalization of Hamiltonian H, Eq. ([3]) 
can be expressed as 

Ak=-E^k.k'^|^^k'- (13) 



This is a self-consistency equation that must be satisfied 
by Ak as a function of wave-vector k and temperature 
T. 



flk 



Uk Vk 



Vv 



u^ 



Qfk 



(4) 



where aj^ and ak are new Bose operators for a single- 
particle excitation of wave-vector k in the superfluid 
state, and coefficients Wk and Wk are so determined as to 
diagonalize the Hamiltonian while maintaining the com- 
mutation rule [ak,Q;jj.,] — <5k,k'- 
The diagonalized Hamiltonian is 



H = Y,iu^. + ^k^ak 



where 



C/k 



^k-^k 



4^k 



(1 + 2nk) ; 



(5) 



(6) 



A. Critical temperature Tc 

Similar to that in the BCS theory of superconductiv- 
ity, energy gap parameter Ak is an important quantity in 
the present theory. It is because of the existence of Ak 
that makes the superfluid state different from the nor- 
mal state. In this and the next subsections we consider 
determination of Ak- 

First, in the limit of T ^ Tc (because of the similar- 
ity between the present theory and the BCS theory of 
superconductivity, we are similarly using T^ instead of 
T\, to denote the critical temperature of the superfluid 
transition), we have |Ak| ^ so that Eq. p^ can be 
linearized and we have an eigenvalue problem: 



Y.y^ 



k,k' 



coth[(ek-Aio)/2fcBTc 



2(ek - Mo) 



(14) 



a 



Ek - M 



(7) 



is the single-particle energy in the normal state, measured 
relative to chemical potential /x; 



^k = da 



(8) 



where /xq is the value of chemical potential /i at crit- 
ical temperature Tc^ and we have used 1 -I- 2rik = 
coth(£;k/2A:BT) and £;k = ek - Mo at T = T^. 

Critical temperature Tc and phase 0k (as defined via 
Ak = |Ak|e*^'') are determined by solving Eq. P^ for 
given interaction Vk.kS single-particle energy spectrum 
Ek and chemical potential /ig. 



Note that it is not necessary to assume Vk,k' < in 
order for Eq. p^ to have a T^ > solution. There- 
fore, the view that an attractive interaction is responsible 
for particle pairing is incorrect. Here we also emphasize 
that the pairing of particles of opposite momenta, as ex- 
pressed by Eq. ([2|), is a kind of ordering in momentum 
space (this point agrees with London's view that super- 
conducting/superfluid state is an ordered state in mo- 
mentum space^ii^); it does not mean that bound pairs of 
particles (due to an attractive interaction) are formed. 



B. |Ak| and £k 

With respect to determination of jAkl, we note that 
the self-consistency equation, Eq. (fT5|) , can be converted 
into 



EiA^ 



d /l + 2nk 



dT 



Eu 



(15) 



by first operating d/dT on Eq. P^ . and then, multi- 
plying the resulting equation by AJ^(1 + 2n\^/2E\^ and 
summing over k. 

Interaction 14;k' no longer appears in Eq. (fT5)) , because 
all information about Vk,k' is already contained in Tc, and 
the latter is involved through the condition |Ak(T'c)| = 0. 

From Eq. (fT5|) we can see that the self-consistency 
equation alone does not allow unique determination of 
|Ak|, because, as one can see, Eq. (|T5|) can have 
an infinite number of solutions. This property of the 
self-consistency equation is true for arbitrary interac- 
tion Vick', because Eq. pS]) is derived for arbitrary 
Vkk'- Actually, this property can also been seen di- 
rectly from Eq. (|13p by noticing that the equation is 
linear with respect to e^^^ . [Even in the case of Vkk' = 
-y, which leads to Ak = ^ Ek'[(l + 2nk')/2^k']Ak', 
|Ak| still cannot be uniquely determined, because phase 
0k' of Ak' in the summation over k' is measured rel- 
ative to phase ^k of Ak on the left-hand side of the 
equation (it is relative phases that matter). This is 
more clearly seen if we re- write the equation as |Ak| = 



^Ek' 



)[(l + 2nk')/2£;k']|Ak'| = l'^Ek'Cos(0k'- 



6'k)[(l^-2nk')/2£:k']|Ak'|, where cos(6lk' - 61k) = 1 or -1 
depending on 0k' — ^k = or tt. We have |Ak| = A = 
k-independent only when 6^ — 6 — k-independent, but 
other solutions for |Ak|, with 0k being k-dependent, are 
also possible. Similarly, for a separable interaction of 
the form Vkk' = ^V^f-^k^k', the solution |Ak| — A|ciJk| 
corresponds to a solution with e*^'' — sgn(a;k)e'^ with 9 
being an arbitrary constant, and is only one of an infinite 
number of possible solutions.] 

On the other hand, we note that diagonalized Hamil- 
tonian H is T-dependent, i.e., both [/k and i?k in Eq. 
^ are functions of temperature T, because of their de- 
pendence upon |Ak(T)|. Since diagonalized Hamiltonian 
H describes a set of independent excitations, and there 
is no transition between different single-particle states in 



thermodynamic equilibrium, we expect the thermal en- 
ergy and entropy associated with a single-particle state 
of wave- vector k to be 



and 



Su 



ek = Uk + UkEk (16) 



-A:B[nklnnk- (l-l-nk)ln(l-l-nk)] , (17) 



respectivelyJ^iiiii^ However, because of the T- 
dependence of H, when we calculate, for a single- 
particle state of wave-vector k, the partition function 

Zk = Trie-"'^/''^'^) (where Hi, = Uk + E^alak), free 
energy _Fk = —/csT In Zk, entropy 5'k = —OFi^/dT and 
thermal energy ek = ^k + TSk, we find that there are 
additional terms involving dUk/dT and dE]i/dT in each 
of the expressions for ek and ^k, as compared to Eqs. 
(dni and ([T7)l . By letting the sum of the additional 
terms in each of the expressions for Sk and S'k to be 
zero, we arrive at 



dT 



"k- 



aEk 

dT 



0. 



(18) 



This equation represents an additional self-consistency 
requirement of the theory, and must be consistent with 
Eq. (1121). To see that this is indeed true, we substitute 
Eq. dll) into Eq. ^ to obtain 



|Aw|2 



df 



l + 2nk 
£^k 



0. 



(19) 



A solution of this equation is certainly also a solution of 
Eq. P^ . and therefore also a solution of Eq. (fTH]) . and 
thus, we see that Eq. (fT8)l is indeed consistent with Eq. 

m- 

Equation (|19p shows that there are two possible solu- 
tions for |Ak| for each single-particle state of wave- vector 
k: one is a trivial solution, |Ak| = (corresponding to 
the normal state), and the other is a non-trivial solution 
(|Ak| > 0, corresponding to the superfluid state) satisfy- 
ing 



_d_ 
df 



l + 2nk 
^k 



which is readily solved to give 
l + 2nk 



E^ 



T-independent for T < Tc- 



(20) 



(21) 



A consequence of Eq. (|20)) . or (|2T|l. is that chemi- 
cal potential /x in the superfluid state is T-independent. 
This is shown in Appendix [X] (where we discuss chemical 
potential /i and number-of-particle distribution (aj^ak) 
in the superfluid state). Then, by using 1 + 2nk — 
coth(£^k/2fcBT) and the condition that |Ak| = at 
T — Tc, we can express Eq. (PT|) as 



(22) 



coth (^k/2fcsT) coth (Ck/2fcBTc) 



El, 



Ck 
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FIG. 1: Energy dependence of energy gap parameter ampli- 
tude |Ak| for —fio/kBTc = 0.5 and different values of T/Tc as 
indicated on the curves. 



FIG. 2: Temperature dependence of normalized energy gap 
parameter amplitude |Ak|/|Ak(0)| for — /.io/^sTc ~ 0.5 and 
different values of ek/fesTc as indicated on the curves. 



with 



a 



Ck — Mo, 



i?k = V(ek-/io)2-|AkP, 



(23) 
(24) 



and fiQ being the value of /i at T = Tc. 

It is shown in Appendix [B] (where we discuss the 
ground state energy of the supcrfluid state) that /io must 
be below a certain negative value in order for the super- 
fluid state to be energeticaUy favorable as compared to 
the normal state. 

Equation (P^ is an implicit solution for |Ak| (or £'k) 
as a function of ek and T for given Tc and ^o ■ A complete 
solution for Ak is therefore a combination of the solution 
of Eq. (122|) for jAkl/fesTc and the solutions of Eq. ((14]) 
for Tc and 0k- 

The present analysis with respect to the determination 
of Ak is similar to that of Ref. !l9 with respect to the 
determination of the energy gap parameter in the BCS 
theory of superconductivity. 

We solve Eq. (122]) numerically by using an iterative 
method^ to obtain lAkl/ZcsTc as a function of e\^/kBTc 
and T/Tf. for given fiQ/ksTc. 

Figure [T] shows ek-dependence of |Ak| for —^o/ksTc = 
0.5 and different values of T as indicated on the curves. 

Figure [21 shows T-dependence of |Ak| for —fio/kBTc = 
0.5 and different values of Ck as indicated on the curves. 

Figure [3] shows Ei^ versus ek for — Mo/^bTc — 0.5 and 
different values of T as indicated on the curves. For 
comparison, normal-state single-particle excitation en- 

(n) 

ergy E^ = ek for T below the Bose-Einstein conden- 
sation temperaturei^iilii^ and E^' = e^ — fiQ at T ^ Tc 
are shown as the dotted curves in the figure. 

From Fig. 3 we can see the existence of an energy 
gap in the superfluid-state single-particle excitation spec- 
trum. Namely, minimum value Emin of _Ek, which is 
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FIG. 3: Superfluid-state single-particle excitation energy E^ 
versus normal-state single-particle energy ek for —^lo/kBTc = 
0.5 and different values of T /T^ as indicated on the curves 
(solid curves). The two dotted curves respectively show 
normal-state single-particle excitation energy -E^ ~ ek for 
T below the Bose-Einstein condensation temperature and 



E, 



(n) 



ek — Mo at r = Tc, as indicated on the curves. 



located at ek = 0, is greater than zero. According to 
Eq. dllD, Era^n = ("Mo) tauh [(-^o)/2fcBTd at T = 0, 
and increases monotonically to Emin = —Mo at T = Tc. 
Since, as shown in Appendix IB| /io must be below a cer- 
tain negative value in the supcrfluid state, we see that 
Emin > in the supcrfluid state. 

Figure m shows T-dependence of E^i^ for different val- 
ues —fio/ksTc as indicated on the curves. 
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FIG. 4: Temperature dependence of minimum superfluid- 
state single-particle excitation energy Emin for —fio/ksTc — 
0.3, 0.5 and 0.8, as indicated on the curves. 



Specific heat C{T) 



FIG. 5: Temperature dependence of specific heat C for 
—fJ,o/kBTc — 0.3, 0.5 and 0.8, as indicated on the curves. 
For clarity, the T/Tc > 1 part of the —fj,o/kBTc = 0.5 curve 
is not shown. The horizontal dotted line shows the value of 
C/iVfcfl = 3/2 for TjT^ > 1. 



Having obtained |Ak(T)|, we can calculated thermo- 
dynamic quantities of the superfluid state. We consider 
the ground state energy of the superfluid in Appendix [B] 
We calculate in this subsection the specific heat of the 
superfluid. 

In the superfluid state, specific heat C is given by 



k 



EllT-^ 



■Z^e^'^/^(l-e-^^/^)2+2(a/r)tanh(Ck/2) 



,(25) 



where we have used Eq. ([22|) . and have adopted a set 
of dimensionless units for the last expression, in which 
energies are measured in units of hsTc and temperature 
in units of T^. 

In the normal state {T/Tc > 1), the specific heat is 
given by 

5nk 



C 



Ea 



dT 



ks 



cCk/T 



^ (e«k/T _ 1)2 

k 



ii 

y2 



^k d{~ii) 
T dT 



(26) 



where we have used the above-mentioned dimensionless 
units for the last expression. 

Figure [5] shows C/Nks versus T/Tc for different values 
of —fio/kBTc as indicated on the curves, where TV is the 
total number of particles of the system, and is given by 



N = 5I"k(T,) 

k 



E' 



1' 



(27) 



where n]f^{Tc) is the Bose function at T = Tc, and we 
have used the above-mentioned dimensionless units for 
the last expression. 

In calculating C/Nks, we have assumed e^ — h k^/2m, 
(where m is particle mass), and have made the substi- 
tution J2]i ^ {2tt)^'^ J (Pk. The integrals involved are 
calculated by using the Simpson method. ?2, The method 
for calculating /i and dfi/dT for the specific heat in the 
normal state is explained in Appendix \X\ 

As shown in Fig. [5l there exists a finite jump in the 
specific heat at the transition temperature, indicating a 
second-order phase transition. The magnitude of the 
jump is larger for a larger value of — /io/fesTc. In the 
limit of T -^ 0, we have dC/dT -^ 0, because of the ex- 
istence of an energy gap in the single-particle excitation 
spectrum. 

Experimentally, the C-versus-T curve shows a A- 
shaped peak at the transitioni^ 



III. 



SUPERFLUIDITY 



We next consider the case where the superfluid is in a 
state of uniform flow with velocity v^ . 

We write the Hamiltonian of the system as 



H = y ' (ck - m) Okflk + 9 E ^kk'a]LaLk«-k'ak', (28) 

"^ kk' 



k 



which is the same as the Hamiltonian of Eq. ([T]) for 
the case of Vs — 0, except that wave-vector k in the 
above expression is now measured in the coordinate frame 
moving with the superfluid. 

We assume that pairing occurs between particles of op- 
posite momenta measured in the coordinate frame mov- 



ing with the superfluid. I.e., we assume 
(a-kflk) 7^ 



(29) 



in the superfluid state for a pair of (k) and (— k) particles. 
Note that, since wave- vector k is measured in the coor- 
dinate frame moving with the superfluid, if we use a free 
Bose gas as an example, a single-particle state of wave- 
vector k means, in the laboratory frame, a single-particle 
state of wave function 



Ok = e 



j(k+q)-: 



and energy 



where 



Ek = ?i (k + q)72m, 



m-Vs/fi. 



(30) 

(31) 
(32) 



Therefore, a pair of (k) and (— k) particles have zero net 
momentum in the frame moving with the superfluid, but 
have a net momentum of 2?iq = 2mwa in the laboratory 
frame. 

Diagonalization of Hamiltonian H is the same as in 
the case of v^. = 0, except that we now have e_k ^ ek for 
Vg ^ 0. The results of the diagonalization are as follows. 

The diagonalized Hamiltonian is 



^-E 



where 



tfk + - (^^k^kak + -B-kaLktt-k 



p(^) 



, (33) 



U\. = -. + -^ + — 73- (1 + "k + n_k) ; 34) 

4 2 4^^^*) 



Ck = sk - m; 



^k -E^ 



is) , ^k - ^- 



(35) 



(36) 



W?^l{^-^±4^-l 



2E, 



is) 



and 



AkUkWk 



2E, 



(«)■ 



(40) 



(41) 



The self-consistency equation for the energy gap pa- 
rameter is 

Ak=-E^>c,>c> '^".";t^" Ak. (42) 



(s) 



2E, 



Since e_k 7^ Ck for v^ ^ 0, we have ^_k 7^ ■Ck, -E_k ^ 



-Ek and n_k 7^ ^k, but we have f7_k 



U^, E^Jl 



E, 



(s) 



and A_k = Ak, as we can see from the expressions shown 
above. 

Chemical potential /i, relative to which energies such 
as ^k and Ek are measured, is q-dependent, because each 
pair of particles in the superfluid state has a net energy 
increase of 2{h q^ /2m) due to the flow of paired particles. 
Namely, we have 



M = A'o + ^ q /2to 



(43) 



where hq is the superfluid-state chemical potential for 
q = 0, and the second term is the per-particle energy 
increase due to the flow of paired particles [see Appendix 
|X]for a detailed derivation of Eq. (|43|) ]. 

For simplicity in presenting the theory, we will use 
the normal-state single-particle energy spectrum of a free 
Bose gas as given by Eq. ([5T|) in the following. 

With the help of Eqs. ^ and (gS]), we have 



and 



where 



^k + ^-k _ ^(0) 
o ^ ^k 



a - e-k h\ 

7, = — kq, 

2 m 



Ao) _ (0) 

?k -*^k "Mo 



(44) 



(45) 



(46) 




is the symmetric part of Ek ; and 



(37) 



and el^ ' is the value of Ck for q = 0. 



|Ak(T,q)| 



nk 



^ LE./knT _ ^ 



(38) 



the Bose function. 

Coefficients Wk and Vk are found to satisfy the following 
relations: 



I ,2 1/ Ck + C-k , , 



(39) 



The following equation is derived as an additional self- 
consistency requirement of the theory: 



1 + "k + JT— k 



e: 



is) 



= independent of T and q. 



(47) 



which is a generalization of Eq. ([21]) to the case of v^, ^ 0. 
We present the details of the derivation of this equation 
in Appendix [Cl 



With the help of Eqs. ([44]) and (gS]) and by using 
l + 2nk = coth{Ei^/2kBT) and the condition that jAk] = 
at (T,q) = (Tc,0), we can express Eq. (|47l) as 



coth 



El' '+2 



2T 



coth 



Ef>-2, 



2T 



2E, 



(s) 



coth 



(0) 



JO) 



A*o 



where 



^r = v(^ 



Moj 



(48) 



(49) 



Zk = cos ak with ak being the angle between k and q, 
and we have used a set of dimensionless units in which 
energies such as E^ , e\^ and /xg are measured in units 
of kgTc, temperature is measured in units of T^, and q 
in units of qq, which is defined via h 0^/2112 = ksTc- 
For ^He, go is ^ 1 A~^, and the corresponding superfluid 
velocity is Wso — ^Qo/"^ ^ 10^ cni/s. 

Equation (|48|) is an implicit solution for |Ak(r, q)| 
or El^'\T,q)]. The variables (k,r,q) for the func- 

in the forms of 



tion |Ak(r, q)| appear in Eq. 

(ej^ ,T, gzk), where gzk is the component of q along k. 

We solve Eq. (|48| by using an iterative method^° to 

obtain |Ak| as a function of ej^ , T and qz^ for given /ig. 

Note that temperature Tck, defined as such that Eq. 

T and qzk if 
cj. and qzk- Only for q = is 
Tck = Tc the same for all single-particle states. 

Similarly, superfluid wave- vector gck, defined as such 



has no |Ak| > solution for given e^ 
T > Tc]i, is a function of e,. 



.(0) 



that Eq. 



(0) 



has no | Ak | > solution for given ej^ , T 

(0) 



and gzk H q > qck, is a function of e^, , T and Zk (or, 
qcWZk is a function of ej^. and T). 

Figure [HJa) shows Tck versus e^, for different values of 

gzk, and Fig. El^b) shows ^ck^k versus e]^ for different 
values of T; —fJ^o/ksTc = 0.5 is assumed in both Fig. 
Elja) and Fig. Efb). 

Note that there is an upper bound for Tck, i-e., Tck < 
Tc, but there is no upper bound for qcu., i-e., qck — > 00 



for 



for ejj. — > 00 , and qck = 00 if k _L q (because Zk 
this case). 

The minimum values of T^k and qck, Tck.miniqzk) 
and qck.,miniT), are of particular importance. For T < 
Tck,min{qzk) [or q < qck,min{T)], the system is in an all- 
paired state, in which |Ak| > for all particles. For 
Tck,min{qzk) <T <Tc [or qck,min{T) < q], the system is 
in a partly-paired state, in which particles in states hav- 
ing Tckiqzk) < T [or qck{T) < q] become de-paired (hav- 
ing |Ak| = 0) while those in states having Tck(gzk) > T 
[or qck.iT) > q] remain paired (having |Ak| > 0). For 







m 



/keT^ 




FIG. 6: (a) Tck versus e^ for different values of qzi^/qo as 
indicated on the curves; (b) Qck^k versus e\^ for different 
values of T/Tc as indicated on the curves; —iio/ksTc = 0.5 in 
both (a) and (b). 



T > Tc, the system is in the normal state, in which 
|Ak| — for all particles. 

A finite viscosity should be observable in a partly- 
paired state, because of the existence of de-paired parti- 
cles, which are expected to behave as normal-state parti- 
cles. We therefore expect critical velocity Vsc{T), defined 
as the superfluid velocity at the onset of an observable 
viscosity, to be about the same as hqck.min(T) / m. From 
the numerical results shown in Fig. [H^b), for example, 
which are obtained for the case of — /io/ftsTc — 0.5, we 
can see that gck,min(T') is a few tenth of go, corresponding 
to a superfluid velocity of a few tenth of Vso = hqo/m. 
Since Vso ~ 10^ cm/s for "^He, we see that the critical 
velocity for ''He in this case is about a few tens of cen- 
timeters per second (which is several orders of magnitude 
smaller than the value previously predicted by Landau^) . 

Numerical results for |Ak| versus ej^ for different val- 
ues of T and qz^ are shown in Fig. [7] [Figs. W[&) and 
[71[b)]. Figure Ela) shows an example of the case of 
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FIG. 7: |Ak| versus e^ for different values of T as indicated 
on the curves; qZk/50 = 0.2 and 0.4 in (a) and (b), respec- 
tively; —^jLo/knTc = 0.5 in both (a) and (b). 



Q < q < (;ck,min(0). In this case, the | Ak]- versus- e 
curve for T = is the same as in the case of q = 
(which is shown in Fig. [T]). As T increases, |Ak 
smaller e^ ' is more strongly suppressed, and decreases 



(0) 

= 
for 



faster. As T increases further so that T > Tck 1 



the 



|Ak|-versus-ek curve has a |Ak| = part for low ener- 
gies (except for e\^' — > 0), for which Tck^qzk) < T [or 
qcw{T) < q]. Figure [7jb) shows an example of the case 
of q > (7ck,min(0). In this case, the |Ak| -versus- e{^ ' curve 
has a |Ak| = part even at T = 0. Namely, at T = 0, 
|Ak| = for those single-particle states with (7ck(0) < q. 

The vertical rises (or drops) in the | Ak|-versus-ek curves 
in both Fig. [Tlja) and Fig. ^h) indicate discontinuities. 
Figures [5] and [5] show, respectively, the numerical re- 
sults for the T-dependence and qzk-dependence of |Ak|. 
As shown in the figures, |Ak| is a monotonic decreasing 
function of T and qz\^, except that, at T = 0, |Ak| > 
is a constant for q < qck(O) (Fig. [5]). The vertical drops 
in some of the curves shown in Figs. [5] and [H indicate 



FIG. 8: |Ak(r, g2k)|/iAk(0, 0)| versus T/T^ for different val- 
ues of e|^ /ksTc as indicated on the curves; qz^/qo = 0.2 and 
0.4 in (a) and (b), respectively; —no/kBTc = 0.5 in both (a) 
and (b). 



discontinuities. 



B. Superfiuid particle density Us 

Particle current density j is the expectation value of 
particle current density operator Jji^ I.e., 



where 



^t Uii] + Uii] * 



the particle field operator 



(50) 



(51) 



(52) 
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and 2.0 in (a) and (b), respectively; —no/kBTc = 0.5 in both 
(a) and (b). 



with 0k (x) being the single-particle wave- function [given 
by Eq. ([50)1 ] , and the velocity operator 



v = -i— V. 
m 



(53) 



A straightforward calculation gives 



j = n— q- -— V(n_k-nk)k, (54) 

m 2m -^ — ' 

k 

where n is the particle density and can be expressed as 
(see Appendix El 



n = y^^nk(T'c), 



n±i^ 



V(4"'-A'o)'-|AkP±^k-q /fcBT 



(55) 



, (56) 



and nk(Tc) is the value of nk at (T, q) — {Tc, 0). 

The first term on the right-hand side of Eq. ^^ repre- 
sents a uniform flow of all particles. The second term rep- 
resents contribution from single-particle excitations and 
de-paired particles, and tends to cancel the first term. 
When all particles are in the superfluid ground state (at 
T — and for q below a threshold), the second term is 
zero. On the other hand, when |Ak| = for all single- 
particle states, the two terms cancel each other so that 

j = o. 

Since superfluid current density j = without pairing, 
as shown above, according to standard quantum theory 
of many-particle systems, it is clear that superfluidity 
arises from pairing of particles, not from Bose-Einstein 
condensation (there is no microscopic theoretical justifi- 
cation for the view that Bose-Einstein condensation leads 
to superfluidity). Therefore, we believe that the super- 
fluid properties of liquid ''He,'^'-'^- as well as the recently 
observed superfluid properties of ultra-cold atomic gases 
(such as the persistent flow of atoms in a toroidal trap^^ 
and the vortices in rotating atomic gasesS^iM), are as- 
sociated with pairing of the atoms involved, not Bose- 
Einstein condensation. 

By using ^k ~ ^ k"^ /2m and making the substitution 
{27r)^'^ J dPk, we can rewrite Eq. ([M]) as 



E 



J = ns—q_ = n^Vs, 
m 



(57) 



where Ug is the effective superfluid particle density, and, 
by using the above-introduced dimensionless units (in 
which energies are measured in units of ksTc, temper- 
ature T in units of Tc and superfluid wave-vector q in 
units of go), can be expressed as 



n 



= 1 - 



ii '^'^°''^"^i '^'"'"^"- 



where 



n 



de 



(0)/^(0)n1/2 







=(E 



(0) 



-Moj 



1 



-nk), (58) 



(59) 



is a function of /io and relates to n (i.e., n ex n), and 
n±k = e^ / — 1 

1. Us for q ^ 

In the limit of q ^ 0, Eq. ([55)1 becomes 

2 r d6^°^(e^''V/'e<°'/T 



(60) 



— = 1- 
n 3nT 



f e^ , 
for q = 0. 



1 



(61) 



where E^ , as a function of e\^ and T for given /iq, is 
determined by Eq. 



1 — I — I — I — I — I — r 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

T/T, 



FIG. 10: Temperature dependence of superfluid particle den- 
sity Tis for the case of g — > for —Ho/kBTc = 0.3, 0.5 and 0.8, 
as indicated on the curves. (The dotted curve shows a linear 
1 - T/Tc dependence.) 



Numerical results for Us/n for the case of g ^- as a 
function of T are shown in Fig. [10] for different values of 
fiQ. Note that dug/dT — > in the limit of T ^ 0, because 
of the existence of an energy gap in the superfluid single- 
particle excitation spectrum. 



2. Us for finite q 

In Appendix \A\ chemical potential /i in the superfluid 
state is determined based on the assumption that all par- 
ticles are paired in the superfluid state. For the case of 
finite q, de-paired particles may exist. The result for /i 
obtained in Appendix|^is no longer valid when de-paired 
particles exist. The question how to determine the chem- 
ical potential in a partly-paired state is not addressed 
in this paper. To proceed, we make the following ap- 
proximation with respect to n±k in calculating n^ when 
de-paired particles exist. For paired particles (for which 
I Ak| > 0), we use the same result for n±k as in the case of 
an all-paired state; and for de-paired particles (for which 
l^kl = 0), we use the following approximation: 



where 






^{n)(rp N| _ ^(eik-A'o) _ I 



n^Tc) 



(62) 
(63) 

(64) 



(in dimensionless units) with e±k — Tt- (±k -I- q)^/2m. 

Numerical results for Us/n obtained based on the 
above-described approximation are shown in Fig. Illf a) 
as Us/n versus T/Tc for difference values of q/qo and 
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FIG. 11: (a) Us/n versus T/Tc for different values of q/qo 
as indicated on the curves; (b) Us/n versus q/qo for different 
values of T/T^, —^o/ksTc = 0.5 in both (a) and (b). 



in Fig. IllT b) as n^/n versus q/qa for difference val- 
ues of T/Tc. In both Fig. [fja) and Fig. [Hl^b), 
^l^o/kBTc = 0.5 is assumed. 



C. Free energy density F 

From diagonalized Hamiltonian H [Eq. ([33|) ]. we de- 
rive the following expression for the free energy density 
in the superfluid state: 



F-^ 



=E 



C/k- 



knT 



■In 1 



2„2 



rq' 

I 

2m 



the 



(e-^ 



first 



Tr e-^/^«^ 



(65) 

term comes from 

which is the usual statis- 



where 
— fc^rin 

tical free energy densityfi^iiiii^ and the second term 
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comes from n(/x — /xq), which is the energy increase due 
to the flow of the superfluid, and which is added to 
F because single-particle energies in the expression for 
Hamiltonian H are measured relative to fi. 

For an isotropic system (such as liquid ^He), F is a 
function of T and q = |q|, i.e., F — F{T,q). As can be 
shown, superfluid current density j and effective super- 
fluid particle density Ug are related to F via the relations 



and 



IdF 
h dq 

m dF 



(66) 



(67) 



respectively [where F is as given by Eq. ([S5)) ]. When 
|Ak| = for all single-particle states, as in the normal 
state, F becomes q- independent, and we have j = and 
n. = 0. 



D. Spatially varying Vs(x) 

The theory presented so far is based on the assumption 
that superfluid velocity v^ is spatially constant. We can 
extend the theory to the case where superfluid velocity Vs 
is spatially varying, by making an assumption that there 
exists a length I such that the following is true: Length 
I is large compared to inter-particle distance so that the 
properties of particles in volume l^ are essentially those 
of an infinite system, but small by macroscopic standards 
so that the volume can be regarded as a "point" macro- 
scopically and all thermodynamic functions of the system 
vary negligibly over the distance I. 

Based on this assumption, quantities such as | Ak|, ris, j 
and F can all be considered as local quantities, obtained 
with respect to particles in a volume l^ around a local 
point X for q = q(x), i.e., |Ak| = |Ak(r, q(x))|, Ug = 
n,{T,q{-K)), j = j(T,q(x)) and F = F{T,q{xj), where 
q(x) = ?7iVs(x)/?i is assumed to vary spatially with a 
length much larger than I. 

The theory presented so far for the case of v^ 7^ is 
similar to the theory presented in Ref. |25| for supercon- 
ductivity in the presence of a magnetic field. 



E. 



Hess-Fairbank effect 



Since ris > in the superfluid state, we see from Eq. 
(1571) that we have 



dF 

dq 



>0, 



(68) 



which shows that a larger value of q is energetically less 
favorable in the superfluid state. This implies that a su- 
perfluid tends to expel superfluid wave- vector q(x), or. 



equivalently, superfluid velocity Vs(x) = ?iq(x)/m, from 
its interior so as to minimize the overall free energy of 
the system. This result qualitatively explains the Hess- 
Fairbank effect,^- the reduction of moment of inertia of a 
rotating cylinder of liquid ^He when it is cooled through 
the superfluid transition. Namely, when the liquid is in 
the superfluid state, because a motionless state is ener- 
getically more favorable, it stops rotating with the con- 
tainer, except in the immediate vicinity of the wall of 
the container where the liquid rotates with the container 
due to interaction between the liquid and the wall at the 
interface. 

Although the theory presented so far provides a qual- 
itative explanation for the Hess-Fairbank effect, it does 
not allow quantitative description of the Hess-Fairbank 
effect. Namely, the theory says that a superfluid tends 
to expel superfluid velocity Vs(x), but it does not tell us 
how Vs(x) can be determined for given boundary con- 
dition and temperature (in the case of a rotating cylin- 
der of liquid ^He, for example, the boundary condition 
is determined by the angular speed and geometry of the 
container) . 

We note that the Hess-Fairbank effect is analogous to 
the Meissner effect in superconductors, and the latter is 
quantitatively describable by combination of the London 
equation^ j = — n^a and the Ampere's law j = V x b 
(here j, Ug, a and b = V x a are, respectively, the elec- 
trical current density, superconducting electron density, 
vector potential and magnetic flux density, and a set of 
dimensionless units is used for the present discussion). 
Our Eq. ([F7|) is analogous to the London equation, with 
Vs playing the role of a. What is missing for a superfluid 
is an equation analogous to the Ampere's law. 

We therefore speculate the existence of the following 
equation: 



where 



j = -AV X uj, 



V X v. 



(69) 



(70) 



is superfluid vorticity, and A a positive constant. 

Note that Eq. ([55)1 applies only to particles (or atoms) 
in the superfluid state; i.e., here j and lj are associated 
with the superfluid component of a fluid, and correspond, 
respectively, to j^ and uJs in a two-fluid model^'^'^'^. 

Equation ([55)1 is clearly only a speculation based on 
the similarity between superfluidity and superconductiv- 
ity, and thus, must derive its validity from experimental 
confirmation of the consequences that it implies. 

Contrary to the common view that superfluid is irrota- 
tional (vorticity-free) (for which there is no microscopic 
theoretical justification), Eq. ([S5)) shows that, analogous 
to that electrical current creates magnetic field, super- 
fluid current creates vorticity. 

As we will see below, an important consequence of Eq. 
(IM)) is the existence of a penetration depth that char- 
acterizes the typical distance to which superfluid veloc- 
ity and superfluid vorticity penetrate into a superfluid. 
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This penetration depth is analogous to the London pen- 
etration depth? that characterizes the typical distance to 
which magnetic vector potential and magnetic field pen- 
etrate into a superconductor. 

We further speculate the existence of an additional 
term in the expression for free energy density F, i.e., 



and 



^-E 



[/^ + :^ln(l-e 



1-e" 



2„2 



+ n 



2m 



Am 



(71) 



where the first two terms are the same as in Eq. (|55|). 
and the third term, which is analogous to the magnetic 
filed energy density in the case of superconductors, is the 
additional term whose existence is speculated. The rea- 
son for this speculation is as follows. We note that Eq. 
([57)1 is an equilibrium property of a superfluid, and thus, 
must also be derivable as a result of the variational prob- 
lem that, in the thermodynamic equilibrium, the overall 
free energy of the superfluid, given by the volume inte- 
gral of free energy density F, is stationary with respect 
to arbitrary variation of v^ (x) . This is true when F is as 
given by Eq. (|7ip. as can be shown with the help of Eq. 

(EH). 

Combination of Eqs. ([57]) and ([S5|) allows quantitative 
description of the Hess-Fairbank effect. In the following 
we present two simple examples. 



1. Superfluid in a rotating cylinder 

We consider a superfluid in a rotating cylinder. As- 
suming the length of the cylinder is much larger than its 
radius R, and neglecting the bottom portion of the cylin- 
der, in terms of cylindrical coordinates (r, cf), z) and unit 
vectors (f , (f>, z), we can write superfluid current density 
j = j{r)cj), superfluid velocity v^ = Vs{r)4> and superfluid 
vorticity u)— uj{r)i, and we have 



Kuj'{r) = ns{r)vs{r) 



and 



v's{r) =^{r) 



Vs{r) 



(72) 



(73) 



where a "prime" indicates a derivative with respect to r, 
ns(r) = ns{vs{r)) is given by Eq. ([58]) . the first equa- 
tion comes from combining Eqs. (j57p and (|69p . and the 
second equation comes from Eq. ([70)1 . These equations 
describe only the behavior of the superfluid component of 
the fluid. We will not consider in this paper the behavior 
of the normal-fluid component of the fluid. 

This is a second-order boundary value problem (which 
is expressed here as a system of two first-order differential 
equations) with the boundary conditions 



v,{R)^Rno 



(75) 



where R is the inner radius of the cylindrical container 
and f2o ttic angular speed of the container. Here we have 
assumed that, in equilibrium, Vs{R) is the same as the 
linear speed of the inner wall of the container, as oth- 
erwise there would be momentum transfer between the 
superfluid and the container. 

In this paper, we will not attempt to solve this bound- 
ary value problem for arbitrary T and Vs{R). Instead, 
for simplicity in presenting the main features of the the- 
ory, we will consider only the case where n^ is spa- 
tially constant. As we can see from Fig. [TlTb). that 
Us is spatially constant is true only at T = for Vs < 
fiQck.miniO) / m [i.e., Us is independent of f ^ at T = for 
Vs < fi.(Zck,min(0)/TO], and is approximately true at higher 
temperatures for sufficiently low values of Vg . 

For a spatially constant n^, the above-described 
boundary value problem can be solved analytically, and 
the solutions are: 



Vs{r) = Rflo 



h{r/\) 
h{R/\) 



and 



w(r) = 



mo /o(r/A) 
A h{R/\y 



(76) 



(77) 



where /„(x) is the modified Bessel function of the first 
kind of order n,^® and 



A = v^ 



(78) 



For R/\ ^ 1, by using the asymptotic expansioi 



,26 



ln{x) ^ e^/v27ra::, where a; ^ 1, we have, near the inner 
wall of the container. 



Vs{r) ~ R^oe. 



-{R-r)/X 



and 



w(r) 



RVto 



-{R-r)/\ 



(79) 



(80) 



from which we see that superfiuid velocity and super- 
fiuid vorticity "penetrate" only a distance of the order of 
A into the superfluid; at a depth of little more than A, su- 
perfluid velocity and superfluid vorticity are practically 
zero; and thus, A has the meaning of "penetration depth" 
that characterizes the distance to which superfluid veloc- 
ity and superfluid vorticity penetrate into a superfluid, 
analogous to the London penetration deptb^ that char- 
acterizes the distance to which magnetic field penetrates 
into a superconductor. 

For R/X ^ 1, by using the approximation^^ Iq{x) ^ 
1 -I- 0(a;^) and /i(x) ~ x/2 + 0(a;^), where x <C 1, we 
have 



w«(0) = 



(74) 



Vs[r) ~ ^oT 



(81) 
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and 



uj{r) ~ 2rJo, 



(82) 



which are the same as the resuhs for a rotating rigid 
body. 



2. Flow of superfluid in a pipe 

We next consider the case where a superfluid flows 
through an inflnitely long pipe of a constant circular cross 
section. Let the axis of the pipe be the 2;-axis, the inner 
radius of the pipe be R, and the total superfluid current 
be /. In terms of cylindrical coordinates (r, ^, z) and unit 
vectors (f , (/>, z), we can write superfluid current density 
j = j{r)z,, superfluid velocity Vs = Ws(r)z and superfluid 
vorticity u}= a;(r)0, and we have 



and 






«sW = -^w> 



(83) 



(84) 



where a "prime" indicates a derivative with respect to r, 
ns{r) = ns(vs{r)) is given by Eq. ([55]) . the first equation 
comes from combining Eqs. ([57]) and (p5)) . and the second 
equation comes from Eq. (j70p . 

This is a second-order boundary value problem (which 
is expressed here as a system of two first-order differential 
equations) with the boundary conditions 



w(0) = 



and 



L^iR) 



I 



2ttRA' 



(85) 



(86) 



■ds. 



where the last condition comes from I = JgJ ■ 

Similar to that for the case of a superfluid in a ro- 
tating cylinder, discussed above, we consider only the 
case where n^ is spatially constant. As mentioned above, 
that Ug is spatially constant is true only at T = for 
Vs < ?i5'ck,min(0)/m, and is approximately true at higher 
temperatures for sufficiently low values of Vg. In this 
case, the above-described boundary value problem can 
be solved analytically, and the solutions are: 



and 



Vs(r) 



uir) 



IX lojr/X) 
2ttRMi{R/X) 

I h{r/X) 
2ttRAIi{R/\)' 



(87) 



(88) 



where /„(a;) is the modified Bessel function of the first 
kind of order n,— and A is as defined by Eq. ([75)) . 



For R/X » 1, by using the asymptotic expansion^^ 
In{x) ^ e^ /\/2ttx, where a; 3> 1, we have, near the inner 
wall of the pipe, 



Vs{r) 



IX 

2-kRK' 



and 



uj{r) ~ - 



-{R~r)/x 



,-{R-r)/\ 



2ttRA 



(89) 



(90) 



from which we see that superfluid flows mainly in the 
region near the wall of the pipe; at a distance of little 
more than A away from the wall, both superfluid velocity 
and superfluid vorticity are practically zero; and A has 
the meaning of "penetration depth" that characterizes 
the distance to which superfluid velocity and superfluid 
vorticity penetrate into a superfluid. 

For R/X <C 1, by using the approximation^^ ^^0(2;) ~ 
1 -I- 0{x^) and Ii{x) ^ x/2 + 0{x^), where x < 1, we 
have 

^s{r) - -4t- (91) 



and 



w(r) 



7ri?2j 



2-kR^A 



(92) 



which show that, in this case, the superfluid flow is nearly 
uniform; and superfluid vorticity is nearly linear in r. 



IV. 



SUMMARY 



We have presented a microscopic theory for superflu- 
idity in an interacting many-particle Bose system (such 
as liquid ^He). The theory shows that, similar to super- 
conductivity in superconductors, superfluidity in a Bose 
system arises from pairing of particles of opposite mo- 
menta. 

In Sec. ini we presented the theory for the case where 
superfluid velocity Vj = 0. The theory shows the exis- 
tence of an energy gap in single-particle excitation spec- 
trum, and the existence of a specific heat jump at the 
transition. 

In Sec. mil we presented the theory for the case where 
superfiuid velocity Vs ^ 0. We derived an equation that 
gives a relation between superfiuid current density j and 
superfiuid velocity v^ (this equation is analogous to the 
London equation for the superconducting state that gives 
a relation between current density of superconducting 
electrons and magnetic vector potential), and an expres- 
sion for superfluid particle density n^ as a function of 
temperature T and superfluid velocity v^. We showed 
that superfluid-state free energy density F is an increas- 
ing function of Vg (i.e., dF/dvg > 0), which indicates 
that a superfluid tends to expel superfluid velocity (i.e., 
a superfluid has a tendency to remain motionless); this 
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result provides a qualitative explanation for the Hess- 
Fairbank effect (which is analogous to the Meissner ef- 
fect in superconductors). We further speculated, based 
on the similarity between superconductivity and super- 
fluidity, the existence of an equation [i.e., Eq. (|69p ] that 
specifies a relation between superfluid current density j 
and superfluid vorticity uj (this equation is analogous to 
the Ampere's law). With the help of this equation, the 
Hess-Fairbank effect can be quantitatively described. 



APPENDIX A: CHEMICAL POTENTIAL ^ AND 
NUMBER-OF-PARTICLE DISTRIBUTION (aj^ak) 

1. yL for Vs = 

We consider in this Appendix chemical potential ^ in 
the superfluid state. We first consider the case where 
superfiuid velocity v^ = in this subsection. 

Chemical potential /i as a function of temperature T 
is so determined such that the number of particles of a 
Bose system is conserved. 

Number of particles TV is given by 



^ = E<4«k) 



which, in the normal state, becomes 



T.[- 



k 



(£k-M)/fcBT _ 



(Al) 



(A2) 



This equation determines chemical potential /x as a func- 
tion of T for given N in the normal state. Here, we have 
assumed volume V = 1 (in arbitrary unit) so that N can 
also be considered as the particle density of the system. 
In the superfluid state, Eq. (|A1|) can be expressed as 



N = 



^l 



(ek - m) 



l + 2nk 



- 1 



(A3) 



by using the results of the canonical transformation of 
Eq. @i and (aj^at) = ^k. 

Since the quantity (1-f 2nk)/-E'k in the above expression 
is T-independent, according to Eq. (HI]), we see that 
particle conservation condition dN/dT = implies that 
dfi/dT = in the superfluid state, which can also be 
expressed as 



/u = /io for T < Tc, 



(A4) 



where fj.Q is the value of the chemical potential at T = 
Tc- Namely, chemical potential fi is T-independent in the 
superfluid state. 

Equation (jA3|) can then be expressed as 



iV = ^nk(T,), 



(A5) 
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FIG. 12: Temperature dependence of chemical potential /i, 
calculated for —fJ-o/ksTc = 0.5, assuming tk = h?k^ /2m. 



For T <n, fJ. 



/io is a constant. The dotted curve shows 
normal state chemical potential fi^"'' for T < Tc\ fi = 
for T < Tc , where Ti ' is the critical temperature 

of Bose-Einstein condensation; and Tc /Tc — 0.49 for 
-fio/ksTc = 0.5. 



where 



ni^iTc, 



A>^k-tia)/kBT^ 



(A6) 



Figure [121 shows the temperature dependence of chem- 
ical potential /i, calculated for —fXo/ksTc — 0.5. In cal- 
culating /i, we have used the Runge-Kutta method. ^° 
Namely, from Eq. (U2|) and the condition dN/dT = 0, 
we drive an expression for dfi/dT, which, together with 
a given value of fi = /iq at T/Tc — 1, allows us to nu- 
merically compute yU (and dpi/dT, which is required in 
calculating the normal-state specific heat in Sec. Ill CI ) 
for arbitrary T by using the Runge-Kutta methodi^ We 
have also assumed e\^ = fi k"^ /2m and made the substitu- 
tion ^^ ex /g°° dee^/"^ in calculating ^(T). The integrals 
involved in the expression for d^/dT are calculated by 
using the Simpson method.^" 

As shown in Fig. [T^l for T > Tc, /i is a decreas- 
ing function of T. For T < 7c, /i = /io is a con- 
stant. The dotted curve in Fig. [T^] shows normal state 
chemical potential a<^"^ for T < Tc] a^^"^ = for T < 
Tc , where Tc is the critical temperature of Bose- 

Einstein condensation;i6,i7,i8 ^^^^ y,jBEC)^j.^ _ ^^^g ^^^ 

— Mo/^bT'c = 0.5. As shown in Appendix [Bl ^q must be 
below a certain negative number in a superfluid state, 
which means Tc > T^?^^, as one can see from Fig. \T% 

Superfluid-state chemical potential /io is an impor- 
tant parameter in the present theory with respect to the 
properties of the superfluid state. For example, zero- 
temperature minimum single-particle excitation energy 
(or energy gap) Emin directly relates to /.to via the rela- 
tion Ernin = (~/^o) tanh[(—/io)/2fcB Tc] according to Eq. 
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((22)) . Generally, /zq relates to the particle density of the 
Bose system, and, as one can see from Eq. (|A5p , a lower 
fiQ (larger |/io|) corresponds to a lower particle density. 



2. (aj^Ok) for v^ = 

The number-of-particle distribution, (aj^ak), is temper- 
ature T-dependent in the normal state. This is no longer 
the case in the superfluid state. The fact that chemical 
potential /i is T-independent in the superfluid state im- 
plies that (aj^Ok) is also T-independent in the superfluid 
state. Namely, from Eqs. (PT|) and (|A4[) we can derive 



d{alak)/dT = for T < T^, 



which can also be expressed as 



{a^au.) = nk(Tc) for T < T^, 



(A7) 



(A8) 



In other words, the number-of-particlc distribution, 
(ajj.ak), becomes frozen at T = Tc when the Bose sys- 
tem is cooled through the superfluid transition. 

Note the difference between the particles described by 
operators a]^ and Ok and the particles described opera- 
tors Qfj^ and Qfk: (akflk) is T-independent (in the super- 
fluid state), whereas (aj^ctk) is T-dependent (for example, 
(aj^ak) ~* asT ~* and {al^au} -^ J^kC^c) as T ^ Tc). 

At T = T^, we have 



(aj^Ok) = (aj^ckk)T=T, = ni^{Tc, 



(A9) 



which shows that all particles in a single-particle state of 
wave-vector k exist as single-particle excitations. 
For < T < Tc, we have 

(aj^flk) = |wkP(akak) + |wkP (l + (a^k'^-k) 

= "k(T,), (AlO) 

which shows that particles in a single-particle state of 
wave-vector k partly exist as single-particle excitations 
and partly are in the superfluid condensate, but the sum 
of the particles remains the same as nk(Tc). 
At T = 0, we have 



(aUk> = kk(0)|' = nk(T,), 



(All) 



which shows that all particles in a single-particle state of 
wave- vector k are in the superfluid condensate. 



3. ^ for Vs 7^ 

We consider in this subsection chemical potential /i in 
the superfluid state for the case where superfluid velocity 



In this case, we have 

^ = E(«k«k) 



jii 



"^k + ^-k / 1 -I- Hk + "-- 



K 



(s) 



- 1 



(A12) 



r(s) 



where ^k and E^^ are given by Eqs. ([35]) and (I37p . re- 
spectively. 

Let X denote any one of T, qi, q2 and q^, where qi 
(i = 1, 2, 3) are components of q = m-Vs/ft- Since in the 

above expression for N the quantity (l+7ik+'T--k)/^k ^^ 
X-independent according to Eq. (I47p . and (^k+C-k)/2 — 
^k + h^q^/2m — /i by using Eq. (I3ip . the particle con- 
servation condition dN/dX — leads to 

d{n^q^/2m - ^i)/dX = 0, (A13) 

which is readily solved to give 

Ai = Mo + fi^q^/2m, (A14) 

where /io is the value of ^l at (T, q) — {Tc, 0), and which 
is the same as Eq. (|43|) . 

Equation (|A12p can then be expressed as 



where 



n'^\Tc) 



^ = E-k (^=) 



g(4'"-/^0)/''-BTc 



(A15) 



(A16) 



and the superscript (0) indicates q = 0. [Equation (jA15|) 
is the same as Eq. ()A5| . as it should.] 



4. (akfflk) for Vs / 

When superfluid velocity v^ j^ 0, as can be shown, 
number-of-particle distribution (aj^ak) in the superfluid 
state becomes 



-,t„. \ _„(o) 



1 



flkflk) = "-k (^c) + 7: (»^k - "--k) 



(A17) 



Ml 



for T < Tek((7Zk) (i.e., for |Ak| > 0), where ^"(Tc) is 
given by Eq. (jAl6p . nk is given by Eq. ([SS]) . and qz^ 
is the component of q = raWg/Ti along wave-vector k (as 
discussed in Sec. IIIIAp . 
We have 



and 



(4ak) + (a^_k«-k) = 2 4°^ (Tc) (AI8) 



(aj^flk) - (ali^a_k) = ?^k - J^-k. (A19) 
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Namely, (flkflk) + (a_k'^-k) is T-independent, but 

(aj^flk) — (aLj^a_k) is T-dependent. 

We have (a]j.ak) = (a_jj.a_k) at T = 0, where all par- 
ticles are paired [for q < (7c,min(0)], and where n\^ — , 



(aj^ak) = (aLkfl-k) = K'k(0)P = nf (T,) (A20) 

for T = 0. 

We have (aj^ak) =/= (a_^a_k) for T > 0, because 

"k 7^ "-k- The difference between (aj^Ok) and (a_j^a_k) 
contributes to the reduction in the effective superfluid 
particle density rig as one can see from the expression for 
ris, Eq. 



APPENDIX B: GROUND STATE ENERGY U{0) 

1. (7(0) for V, = 

The ground state energy of the superfluid state in the 
case where superfluid velocity Vs = is 



C/(0) 



(H) 



T=0 



E^i'W 



-E' 



a 



Mi^/ksT, _ y 



(Bl) 



where ^k = Ck — Mo and we have used -Ek(O) = 
fktanh(^k/2fcBT'c), according to Eq. (g^)- 

It is interesting to note that the right-hand-side of Eq. 
(|B1I) can be interpreted as the thermal energy of a Bose 
system of paired particles at T = T^ having a pair exci- 
tation spectrum of 25k- 

Since the single-particle energy in the expression for 
Hamiltonian H is measured relative to chemical potential 
/I = /zo, we must add the term /zo-/V to the above expres- 
sion for C/(0) when we compare C/(0) with [/(^^^) = 0, 
the energy of a Bose-Einstein condensate, which is the 
ground state of the normal statefi^i^ii^ Namely, in order 
to have a superfluid state, we must have 



c/(o) + ^ioN 



E 

k 



6 



Mo 



s2Ck/feBT, _ I git^/kBT, _ I 



< 0, 



(B2) 



which shows that /ig must be below a certain negative 
value /i* in a superfluid state. Assuming ek = fi- k^ /2m 
and making the substitution J^k ~^ {1/2tt)^ J (Pk, we 
find [from the condition that U{Q) + fiQN = at /iq = /i*] 
that ti*/kBTc ~ -0.21. 



2. f/(0) for v, / 

When superfluid velocity v^ ^ 0, the ground state en- 
ergy is 

C/(0) = C/(o) (0) + {h\^/2m)N, (B3) 



where t/^°^(0) is the ground state energy for Vs = 0, 
as given by Eq. (|Bip . and the second term comes from 
(/i — Mo)A^, which is the energy increase due to the flow 
of the superfluid (here /i is the chemical potential in the 
superfluid state when Vg ^ 0, /iq is the chemical potential 
in the superfluid state when Vs — 0, and h q^/2m = 
mv'^/2 is the per-particle energy increase due to the flow 
of the superfluid). 



APPENDIX C: DERIVATION OF EQ. (|47|) 

We present in this Appendix the details of the deriva- 
tion of Eq. (|47|. 

For convenience, we define 



Ck = 



1 + JT-k + n-k 



2E, 



is) 



Then, the self-consistency equation, Eq. 
rewritten as 



Ak = -E^kk'Ck'Ai, 



(CI) 



can be 



(C2) 



When superfluid velocity v^ ^ 0, we expect Ak to be 
a function of temperature T and superfluid wave-vector 
q = m'Vs/h. Let X denote any one of T, qi, q2 and 93, 
where qi {i — 1,2,3) are components of q. We operate 
d/dX on both sides of Eq. (jC2|l to obtain 



aAk 

dX 



E^^ 



dCk, dAk, 

k' ^,^ Ak' + L-k' 



dX 



dX 



(C3) 



(note that Vk,k' is assumed to be X-independent). 

We next multiply both sides of the above equation by 
CkAJ^, and then take summation over k, i.e.. 



VCkAJ^ 



k 



9Ak 
dX 




K^^""*^ 



(C4) 



The quantity inside the first pair of parentheses on 
the right-hand side of the above equation equals to A^, , 
according to Eq. (|C2p . Therefore, the second of the two 
terms on the right-hand side is the same as the term on 
the left-hand side. Thus, we have 



EiA. 



.9Ck 

dX 







We want a jAkj > solution. Clearly, 
Ck = A-independent, 
which is Eq. (f47|) and satisfies 
5Ck 



dX 



0, 



(C5) 



(C6) 



(C7) 



17 



is a solution of Eq. (|C5p . 

However, Eq. (|C6[) is not the only possible solution of 
Eq. (|C5P [as one can see, Eq. (jC5j) actually can have 
an infinite number of solutions]. We therefore need to 
justify that Eq. (|C6p is the only physical solution, which 
is done as follows. 

Diagonalized Hamiltonian H of Eq. ([55)1 describes a 
system of independent quasi-particle excitations. In ther- 
modynamic equilibrium, there is no transition between 
different quasi-particle states, except pairing correlation. 
Therefore, we can calculate, for each pair of (k, — k) ex- 
citations, the partition function Zj. = Tr(e^^''/'^^-^), 
where H\^ = 2Uk + E^al^ak + E-j^^a ^^a^i^, free energy 
Fk — —A:bT In Zk, entropy 5k = —dF\s./dT and thermal 
energy Ek = ^k + TS^^. 

We expect entropy S'k and thermal energy ek to be 

S'k = -fcs [?iklnnk - (1 -l-nk)ln(l + n-k) 

+ n_klnn_k- (l + '^-k)ln(l+n_k)] (C8) 



and 



£k = 2C/k + "k^'k + Ji-k^'-k, 



(C9) 



respectively. 

is) 

However, because C/k and El. ' are T-dependent, there 
are additional terms involving dU\^/dT and dE^' /dT in 
the expressions for S'k and £k obtained from i^k , as com- 
pared to Eqs. (jCSp and (|C9|) . By letting the sum of the 
additional terms in each of the expressions for S'k and ek 
to be zero, we arrive at 



. 9t/k 
^ dT 



(^) 



("k 






(CIO) 



Similarly, we expect the contribution from each pair of 
(k, — k) excitations to the particle current density of the 
superfluid to be 



jk = {n-k - n_k) — k 
m 



(Cll) 



(where fik/m is the velocity of an excitation in the state 
of wave- vector k) . 

We further expect jk — h~^dFk/dq. However, as 
compared to Eq. (jClip . the expression for jk obtained 
from Fk contains additional terms involving dUk/dq and 

dE^' /dq. By letting the sum of the additional terms in 
the expression for jk to be zero, we arrive at 



is) 



^dUk , , dE, 



= 0. 



(C12) 



Substituting the expressions for U^ and E^^' [Eqs. 
and dSll)] into Eqs. (jUIO)) and ((CT2)) . we obtain 



lAkl^^^O, 



(C13) 



and therefore, for |Ak| > 0, Eq. jCll [which is Eq. (|47)) ]. 
The derivation of Eq. (j47p presented in this Appendix 
is similar to that presented in Ref. '25! for a similar equa- 
tion in the theory for the superconductivity in the pres- 
ence of a magnetic field. 
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